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Self-Consistent Theory of the Gain Linewidth for Quantum-Cascade Lasers 
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The linewidth in intersubband transitions can be significantly reduced below the sum of the 
lifetime broadening for the involved states, if the scattering environment is similar for both states. 
This is studied within a nonequilibrium Green function approach here. We find that the effect is of 
particular relevance for a recent, relatively low doped, THz quantum-cascade laser. 



Material gain engineering in semiconductor quantum 
wells is a major topic in the physics of semiconductor 
laser emission [T|. In particular, the position of the peak 
gain and the linewidth of the active medium are of im- 
portance for applications in telecommunication or short 
pulse generation. Whereas the gain properties of semi- 
conductor interband quantum well lasers are well un- 
derstood 0], intersubband emitters, such as quantum- 
cascade lasers (QCLs) |3J, are still a subject of inten- 
sive research. Major differences of QCLs in comparison 
to interband lasers result from the similar in-plane band 
curvature for the electronic transitions and the low en- 
ergy collective excitations involved in the photon emis- 
sion. Besides model calculations 0, 0| for a standard 
QCL design a detailed many-particle study of the 
gain linewidth of intersubband lasers for the broad vari- 
ety of QCL structures available 0, M, E| is still missing. 
The purpose of this Letter is to evaluate a self-consistent 
theory for the influence of the electronic scattering mech- 
anism on the gain linewidth of QCLs and to clarify its 
impact on the description of QCL structures ranging from 
the infrared to THz-regime. 

For interband lasers the many-particle interactions are 
well described on the level of the second order Born ap- 
proximation Due to the expansion of the scattering 
matrix in wavenumber states, the linewidth of the tran- 
sitions is determined by diagonal (proportional to the in- 
verse lifetime) and non-diagonal dephasing mechanisms 
[Tol ITl) . The consistent description of the detailed inter- 
play of diagonal and non-diagonal scattering mechanisms 
[l2j results in a proper prediction of the gain shape. In 
particular spurious absorption below the band gap could 
be eliminated from the theory |l3j . Our aim is to evalu- 
ate similar information for the gain in the intersubband 
gain of QCL structures. So far, only simple two band 
intersubband emitters 0, 0] have been investigated on 
this self-consistent level. 

Here, we bridge the gap to real QCL devices, and ex- 
plore the influence of a self-consistent description of the 
electron-phonon and electron-impurity scattering on the 
gain linewidth. Although excluded in the present study, 



electron-electron scattering can be treated on the same 
footing 01 . Our results are illustrated for two different 
samples: (A) the terahertz QCL reported in Ref. (B) 
a midinfrared QCL reported in Ref. [||. We demonstrate, 
that the different behavior can be understood by the im- 
pact of non-diagonal dephasing on the linewidth. Finally, 
we give general conditions, when non-diagonal dephasing 
is important and thus a full quantum theory of the gain 
is indispensable. 

First, we evaluate the nonequilibrium electron occu- 
pations with nonequilibrium Green functions Q using 
nominal sample parameters. Here, however, we extend 
the model in two ways: (i) Scattering at ionized impuri- 
ties is treated microscopically using Debye screening by 
a fictitious homogeneous 3D electron gas matching the 
average doping density, (ii) We include the non-diagonal 
elements in the self-energies, so that Eqs. (6-8) of Ref. 
take the form: 
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and similarly for the electron-phonon scattering processes 
(A denotes the normalization area and m c is the effec- 
tive mass). The values of 7cuct 2/ 3 1/ 3 2 are taken for a typical 
fixed momenta |k| , |k' | , while the integration over Z (k, k') 
is performed. The current voltage characteristic is eval- 
uated along the lines of Ref. [j an d is shown in the inset 
of Fig. □ for the THz-QCL A|16|. It agrees excellently 
with the data presented in Fig. 2 of Ref. . 

Next, we evaluate the gain according to Ref. [13. The 
material gain coefficient g(fl) reads as a function of the 
frequency f2 |2|: 
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where e r = 13 is the dielectric constant. The change in 
the current density 5 J for a perturbative electric field 
SF(t) = SF{n)e~ int is given 



by 0: 
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Figure 1: Calculated gain at an operating bias of 11.5 V for 
sample A (Ref. The inset shows the current- voltage char- 
acteristic where the bias equals 177 times the bias drop per pe- 
riod Fd and the current is evaluated for an area of 54000/im 2 . 
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Figure 2: Calculated gain for sample B (Ref. 

Neglecting the (5£-terms provides us with, see Ref. |l7) : 



where the states ^i(z) are the eigenstates of the 
heterostructure potential including the self-consistent 
Hartree field, Zij is the dipole matrix element, and d is the 
period of the structure. The field-induced changes 5Gf 3 - 
of the Green functions are evaluated self-consistently in- 
cluding the non-diagonal elements of the self-energies 
in the terms (full theory). The results are given in 
Figs. 11121 (crosses) for the two samples considered here. 
For comparison, we have evaluated the gain by the simple 
approach (using lifetime broadening only): 
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Here, we evaluate the level energies Ei and their re- 
spective widths Ti (FWHM) from the spectral functions 
— 23{G^ ct (fc = 0, E)} and the densities are given by 

m = £fc/d^{G«(fc>#)}/M)- Eq- @ can be de- 
rived restricting to diagonal dephasing, see, e.g., Ref. [l8| . 
where effects due to different effective masses were stud- 
ied, which are neglected here. It has also been used in 
Ref. Q except for the replacement Ei — Ej = Kl in the 
first line and the neglect of counter-rotating terms with 
Ei < Ej. Figs. 1 1121 show, that the simple approach is in 
good agreement with the full theory for sample B, but 
exhibits a far too large width of the gain peak for sample 
A. This shows, that the simple (Lifetime) approach docs 
not correctly describe the gain spectrum for all QCLs. 

In the following we want to shed light into the physi- 
cal mechanisms which cause the failure of the simple ap- 
proach For this purpose we assume that the Green 
functions Gij of the stationary state are diagonal in our 
basis of eigenstates and given by their generic forms 
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<5G<°(k, E) = G?T(k, E + nn)(-eSF{n)z i:j )Gf(k, E) 
+ G<(k, E + htyi-eSFityz^Gf v (k, E) (7) 

We define Sg<"(E) = f™ dE k SG^°(k, E). Extending 
the Ek integration to — oo (which is possible only for 
E > Max{£, - Utt, Ej}) we find 

^ {E) ~ 2m ^ + M7-£ 4 + i(r J + r,)/2 (8) 

which justifies Eq. 10} by inserting into Eqs. IL'l.'it . Thus 
the simple approach can be related to the neglect of 5E- 
terms in our full theory. 

Now we want to examine the ST, < -terms (neglecting 
^rct/adv £ Qr p Ur p 0se f simplicity in the analytic 
approach), which give: 

5G<{k, E)^5G<°(k, E) 

+ Gf (k, E + hSl)6X<(E)Gf v (k, E) (9) 

For i 7^ j, the dominant contribution [2(j to <5£^' lmp 
takes the form 

/>00 

«5£<' imp (£) = 7 ™ p /2^ / dE k ,5G<(k>, E) , (10) 
Jo 

which gives the non-diagonal depha sing as a k' integral 
similar to density matrix theory \wi fllL Ha . Il3| . Now we 
define Sg^(E) — J dEfc5Gy(k, E). Again extending the 
E k integration to — oo, Eqs. I|5I1U|) give 
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which has the form discussed in Ref. [19|. We immedi- 
ately see, that the width of the gain peak is reduced. 
This relates to the line narrowing due to non-diagonal 
elements in the self-energies (or non-diagonal dephasing, 
see also 0,^^)- Impurity scattering contributes with 
approximately 7™ p to the full broadening of the states 

and 



i. Therefore it is crucial to compare 7„'"; P with 7 
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Figure 3: Band diagrams for samples A and B together with 
main levels. 
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From Eq. Q we find that 7S? is of the order of 



luiiljjjj if the scattering potentials V™ v and V-™ p 
are correlated, i.e. both states ^i(z) and *&j(z) are ex- 
posed to an identical scattering environment. In this case 
the transitions are less affected by the scattering, result- 
ing in a narrowing of the gain feature. Indeed we find for 
the dominating 4-5 transition (see Fig. |3J) for sample A 

7^4 P 5 = 2-2 meV, 7^44 = 2-7 meV, and 7^5 =2.2 meV. 
This means both scattering environments are highly cor- 
related, resulting in the strong narrowing in the gain 
visible in Fig. ^ In contrast the main gain peak at 
eFd = 260mcV in sample B can be attributed to the 
7-8 transition with 77™^ = 1.1 meV, 77777 = 1.9 meV, 



and 7g™8 5 8 = 8.3 meV. Thus, correlations in the scatter- 
ing environment are weak and the simple approach works 
well, see Fig. 

The difference between the two designs can be under- 
stood by considering the nature of the scattering poten- 
tials. Sample A is low-doped, corresponding to a screen- 
ing length of 30 nm, and the differences between the 
wavefunctions 4 and 5 occur on a smaller scale, see Fig.|3 
In contrast the screening length is 8 nm for the high- 
doped sample B. Thus the scattering potential is much 
more local for sample B causing significant differences in 
the scattering environment for both levels. This effect is 
strengthened by the fact that interface roughness (with 
5-function potentials) is also of relevance for sample B. 

Note that the analytical motivation given above repro- 
duces the right trend, but cannot be used for quantitative 
analysis as: (i) Phonon scattering shows also some line 
narrowing due to correlations in the potentials, which 
is contained in the full theory, (ii) The influence of 
<5£^ t/,adv (E) has been neglected, (iii) Significant non- 
diagonal elements in Gfj (E) are present even in the basis 
of energy eigenstates. 

In conclusion we have shown that the simple approach 
|@J for QCL gain is not reliable if the states involved 
in the lasing transition are exposed to the same scatter- 
ing environment. Non-diagonal dephasing then becomes 
important and narrows the linewidth below the sum of 
the individual widths of the levels (lifetime broadening). 
This requires a full quantum kinetic description of the 
gain in QCLs. 
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